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Abstract. We present here the solution of the problem on linearization of fourth- 
order equations by means of point transformations. We show that all fourth-order 
equations that are linearizable by point transformations are contained in the class of 
equations which is linear in the third-order derivative. We provide the linearization 
test and describe the procedure for obtaining the linearizing transformations as well 
as the linearized equation. 
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1. Introduction 

The problem on linearization of second-order ordinary differential equations by means 
of point transformations was solved by Sophus Lie |T] in 1883. More specifically, he 
showed that the linearizable equations are at most cubic in the first-order derivative 
and gave the linearization test in terms of the coefficients of these equations. 

In 1997, G. Grebot [2] studied the linearization of third-order equations by means 
of a restricted class of point transformations, namely t = ip (x) , u = ip {x, y) . However, 
the problem was not completely solved. 

In 2004, N.H. Ibragimov and S.V. Meleshko \B\ solved the problem of linearization 
of third-order equations by means of point transformations. They showed that all third- 
order equations that are linearizable by point transformations are contained either in 
the class of equations which is linear in the second-order derivative, or in the class 
of equations which is quadratic in the second-order derivative. They provided the 
linearization test for each of these classes and describe the procedure for obtaining 
the linearizing transformations as well as the linearized equation. 

The present paper is devoted to obtain criteria for a fourth-order equation to be 
linearizable by change of the dependent and independent variables. In our calculations 
we used computer algebra packages. The final results were checked by comparing with 
theoretical results on invariants as well as by applying to numerous known and new 
examples of linearization. The paper is organized as follows. 

2. Point transformations of fourth-order equations 

We consider the fourth-order ordinary differential equation 



to equation ([T]). 

We begin with investigating the necessary conditions for linearization. The general 
form of ([T]) that can be obtained from linear equations by any point transformations ([2]) 
is found on this step. In consequence, we identify two candidates for linearization. 

A linear fourth-order ordinary differential equation we use in the Laguerre form. In 
1879, E. Laguerre showed that in linear ordinary differential equation of order n > 3 the 
two terms of orders next below the highest can be simultaneously removed by equivalence 
transformation (see [4], Section 10.2.1 and the references therein). Therefore, we write 
the general linear fourth-order equation in Laguerre's form 




(1) 



We apply a point transformation 




(2) 



u 



+ a{t)u' + I3{t)u = 



(3) 



where t and u are the independent and dependent variables, respectively. 
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2.1. The candidates for linearization 

Considering t and u as the new independent and dependent variables, respectively, one 
obtains the following transformation of the first-order derivative 

, _ D^{ip) _ ipx + y'ipy 



where (Px = ^y = ^^c., and 



(4) 



d 



d 



D^ = — + y'— + y"— + y 



d 



+ y 



(4). 



d 



+ 



dx dy dy' dy" dy'' 
is the total derivative. Likewise, one obtains the transformation of derivatives of the 
second and higher order. Namely, denoting by P{x,y,y') the right-hand side of (jl]), 

i'x + y'ipy 



P{x,y,y' 



one has 



u 



(fx + y'(py 



A 



Dx iv) ^x + y'^y {(fix + y'vyY 

Denoting by Q{x,y,y', y") the right-hand side of ([5]), 

A 



Q 



{^x+y'^y)' 



y" + 



one has 



u 



A 



x+y'Qy+y"Qyt+y"'Qy„ 

'Px+y'ipy 

{^x + y'^y)y"' -^^y{y"f 



{fx+y'fyY^ 

Denoting by y, y' , y", y'") the right-hand side of 

A 



R 



{if^ + y'ipyY 



{fx + y'fy) y'" - ^fy {y'f 



+ 



+ 



hence. 



u 



(4) 



D^Cfl) ^ Rx+y'Ry+y"Ry, +y"'Ry" Ry'" 



(ip^+y'ipyY 



+ 



Thus, (El) becomes 



{fx + fyy'fAy^"^ + [-10A(¥;, + fyy')fyy'' 



{fx + y'fyV 

-2(2(5v3a;yA - iPyA^)iPy + {^fyyA - Aif yAy)(p ^)y''^ 

-"^{^{^fxyfx + fxxfy)^ " "^{fx^^y + '^fy^x)fx)y' 



+ • • -v + 



0. 



Here 



(5) 



(6) 



(7) 



(8) 



A = fx'^y - fy-lpx 7^ 
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is the Jacobian of the change of variables ([2D- It is manifest from ([H]) that the 
transformations ([2]) with yj^ = and (py ^ 0, respectively, provide two distinctly different 
candidates for linearization. 

\i iPy = Q we work out the missing terms in ([8]), substitute the resulting expression 
in ([3]) and obtain the following equation 

^(4) ^^AW + A,)y"' + B,y"' + {C^y'^ + C,y' + C,)y" 
+D^y'^ + D^y'^ + D^y'^ + Diy' + = 0, 

where 



(9) 



Aq = -2{Lp^1Py)-^{3Lp^^iPy - 2Lp^i)^y), 
Bo = S{^/Jyy^tjjyy, 
C2 = 6{i)y)~'^i'yyy, 







Co = 






D,= 




D3 = 


-2{^^i)y) 


(J^'^xx^yyy 


D2 = 




[^^xxx^x^l 


Di = 




^ipvlxi^y 



■'xyyj} 



Q{?>'Pxx'^xy - 'fx^xxy)^x 



^xyyy)-, 
„2 „/, 



^xxyy ) 5 



^xxy I 



^xxxy j 



-2(hLp^^1py - 4:(p^i!^y)ip^^^(p^ + IfixxxxV^l-ipy 



(10) 

(11) 

(12) 
(13) 
(14) 
(15) 
(16) 
(17) 
(18) 

(19) 



Do = -{^l-ipy 



(ISV^L - V^x« + ^xxxx^l)'ipx - {W^xxxV^xx-ipx - ^^xxx^x^x 
+ 15v?LV'xz - Q^xx^xi^xxx + ^ll3i' + ^llpxxxx)^x ■ 

Definition 1. We call ([9]) with arbitrary coefficients Aq = 
Ai = Ai{x,y), Bo = Bo{x,y), Co = Coix,y), Ci = Ci{x,y), C2 
and Di = Di{x,y), {i = 0, . . . ,4), the first candidate for linearization. 

If ipy 7^ 0, we proceed likewise and setting r{x,y) = arrive at the following 
equation 



(20) 

Aix,y), 
C2ix,y), 



y^'^ + YTr^-W + F^y'"" + Fry' + Fo)y" 
[Ihy"^ + {H^y'^ + H^y' + Ho)y"^ 



+ 



+W + Jzy'^ + J2y'^ + Jiy' + Jo)y" 

+Kry" + K,y'^ + K,y'' + K^y'^ 
+Ksy'^ + K^y'^ + K^y' + Ko] = 0, 

where 

F2 = -2{iPy/\)-\biPyy/\ " 2iPy/\y), 



(21) 



(22) 
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Fi = 4((^j,A)-i [(A^ + A^r - 5rj,A)(^^ - hipyyr/S\ , 

Fq = -2{ipyA)-' [{{5ryA - 2A^)r + 5r^A)^y + 5^yyr^A\, 

H2 = QiiPyAy^idifyyA " 2iPyAy), 

H, = -3{ipyA)-^ [(5A^ + 3Ayr - 25ryA)cpy - 2Q(pyyrA\, 

Ho = Si^yA)-' [(5(3r^ + 2ryr)A - (5A, - Ayr)r)^y + lO^yyr^A] , 

J4 = -{iplA)-\lOiPyyyiPyA " 45(^jj,A + SOlfyylfyAy - GlflAyy), 

J3 = 2(^jA)-i [3((2(A., + AyyT - 5ryAy) - 5ryyA)^l 

-5((A^ + 3Ayr - iryA)(Py - 6(PyyrA)iPyy) " 20(p yyylf yV A] , 

J2 = G{iplA)-^ [(A^^ + Ayyv^ + 4A^yr - 5(2A^ + 3Ayr - 5ryA)ry 

-lOvyyvA - 5r,Ay - 5r,yA)^l - 5(((3(A, + Ayr) - 10ryA)r 
-2r^A)ipy - 9(pyyr^A)ipyy " lO^^j/yj/V^j^r^ A] , 

Ji = -2((^J^)-i [((5(3(3A, + Ayr) - 14r,A)r, - 6(A,,r + A,,) 

+20ryyrA)r + 5(3(A^ + Ayr) - 16ryA)r^ + 5r^^A + 20r^yrA)ipl 
+15(((3Aa; + Ayr - 8ryA)r - 4rxA)(py - 6ipyyr'^A)ifyyr 

+ 20iPyyyiPyr^A , 

Jo = -{<pIA)-^ [((2((5r^^rA - 3A^^)r + 5r^^A + Sr^^rA) 

-5{7ryA - 6A^)ryr)r - 5{2{7ryA - 3A^)r + 9r^A)r^)ipl 
-5{3{2{{2ryA - A^)r + 2r^A)ify + 3^yyr^A)^yy 

-2lfyyylfyr'^A)r-' 

^7 = -i^flA)'^ [Pyyyyfli^y - ^^VyyyVyyVyi^y + ^Vyyy^l^^yy + ^^^lyi^y 
-i^fly^'yi'yy + Q¥>yy¥>lAyy " " V'JV'y" " ' 

^6 = (V'yA)-^ [3(5((7(/?yV2/2/^ - ^\)Vy - ^{(Py^yr - A)(pyy)(pyy 

-2{7(Pylljyyyr " 5 Ayy) (fl) (f yy + {llfilfji/jr + llfi^i/jyar " (/J^ftA 

+ 7(/?2/V'2/2/ra'^ - 4:Ayyy)(fil + 2{35(fi yyifi y^^yr " 30(fiyyA " UifllPyyr 

+ 10iPyAy)iPyyyiPy - yi^ yr " 5 A) (^yy^/^ ] ' 

^5 = [(2(3(A^j,2, + 3Aj/2,j^r - hryAyy - hryyAy) - hryyyA) 

-3(7(/?J/3V'r + l^li^yar - 2(plaA + 7iljyyyyr)ipyr)ipl 
-3(2(5(A^2/ + ^\yr - ^T'y^ - '^^yy^) - '^'^'Pyi^yyy^Vy 
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-15((A3. + llAj,r - 3rj,A - 7(pyijjyyr'^)(py 

+ 7{iPyljjyr - 2A)iPyyr)iPyy)iPyy " 2((5(A,; + H AyT " SVyA) 
-A2iPylljyyr^)iPy + 15(7^2^^2/'^ " 1 2 A) (^J/^/'^) (^^/TO^J/ 

+ 3{7iPylljyr - 10A)(Pyyyyiplr] , (SS) 

K4 = -((^jA)-i [(2(45rj,j,rj,A - lOvyyA^ - SSr^/^^Aj^r + 50rjAj, 

2iQryAxy 50TyAyyr -\- iiA^yyT "I" QiAxxy ~\~ ^'^^yyy^ 

^OVyyyVA ^TxAyy XOVxyAy ST^jj^j^A) 

-5{7iplPipr + 7cplipyar - 3(plaA + 7'ipyyyyr)(fiyr^)(fil 
+15((3((5(A^ + 5Aj,r) - UryA)r - r^A) - 35ipyjjjyyr^)ipy 
+35{(pyijjyr - 3A)(pyyr^)(ply - lO{Axx + 3lAj,2/r^ + 13Aa;j,r 
-8(A^ + 6Ayr - 2ryA)ry - 26ryyrA - AvxAy - AvxyA 

-2liPyll)yyyr^)iPyyipl " 1 ( ( ( 5 ( A^ + h AyV) " UT y A)r " A 
-UiPyi^yyr^)iPy + 5(7(^j;'0j;r " A) if yyV"^) yyy^f y 

+ 5{7(Pylljyr - 15A)(Pyyyy(plr^\ , (36) 

K3 = -i^lA)-^ [((13A^^j, + 35Ayyyry + Axxcc + ^l^^yyV^ 

-5{3Axx + 26Ayyr'^ + 23Axyr - (15A^ + 49Aj,r - 25ryA)ry)ry 
— 5(13Aj; + 32Aj^r — 50ryA)ryyr — 65ryyyr^A — 5(3Aa;j, + 5Ayyr 

\QTyAy 7TyyA^Tx ^"^xx^y ^"^xxy^ 

-5(3Ax + llAyT - 15ryA)rxy - SOvxyyvA - 5(7ippijjr + 7ipl'^yar 
-AiflaA + 7i^yyyyr)ipyr^)ipl - 5(2((2(2A,, + ITA^^^r^ + llA,^r) 
-(29A^ + 75Ayr - 51ryA)ry - A5ryyrA)r - (3A^ + 13Aj^r - 13ryA)r^ 
-rxxA - UrxyrA - 21^ yi,yyyr^)^l - 3((6((5(A^ + 3A^r) - 13r^A)r 
-2r^A) - 35ipyijjyyr^)ipy + 35(v?j/^j,r - AA)^yyryyyr)(fyy 
-10(2((5(A^ + 3Aj,r) - 13ryA)r - 2r^A - 7^y^yyr^)^y 

+ 5{7iPylljyr - 24:A)(Pyyr^)(Pyyy(Pyr + ^{7^y'4^yr " 2Q A)^ yyyy^lf^^ , (37) 

K2 = -i^lA)-^ [((3((5A^^j, + 7Ayyyry + Axxx + 7Axyyr^) 

-(3(13A^^ + 28A^j,r' + 39A^^r) + (204rj,A - 161A^ - 217 Ayr)ry)ry 
-{79Ax + lldAyr - 26AryA)ryyr - 5Aryyyr'^A)r 
-{S{2Axx + 7Ayyr^ + llAxyv) + {171ryA - 64A^ - UOAyr)ry 
-72ryyrA - 18rxAy)rx - (4A^ + llA^r - 21ryA)rxx - l2rxxyrA 
-Txxx^ - {{37 Ax + 53Aj/r - 150rj^A)r - 33rx.A)r:,;.y - SSr^^^^r^A 
-3{7(pl(3^r + 7^li)yar - SyjJaA + 7i)yyyyr)^yr'^y^l 
-3(2(5((2A,, + 7Ayyr^ + 6A,,,r - (13A, + 19A,,r - 20r,A)r^ 
-13ryyrA)r'^ - {{3Ax + hAyV - llryA)r - rxA)rx - TxxtA 
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-Qr^yr^/\) - 2lipyijyyyr'')ipl " 15((2((5(A, + 2Aj,r) - 12r^A)r 

-3r^.A) - 7^yll)yyr^)^y + ^ y^) yT " ^ /\) ^ yyV"^) ^ yyV"^) ^ yy 

-2(2(5((5(A, + 2Aj,r) - l2ryA)r - 3r,A) - 2l^yijyyr^)^y 

+ 15{7ify1Pyr - 30A)lfyyr^)(fyyy(fyr^ + 3 {7 (f y-^ yV " 2 5 A ) (/^yj^yj^ ^5 , (38) 



Ki = -iflA)-^ lii^i^xxy + Ayyyry + 3 A,,, + 7A, 



xyy' 



(33A^^ + 28Ayyr^ + 49A^j^r + 2{59ryA - 56A^. - 42Ayr)ry)ry 



-(43A^ + 42Aj^r - 128ryA)ryyr 



r^A)r^ 



21Ayr)r^)r^ 



9^xxy^ A 



'^f xxx'^ A 



-((12A^^ + lAyyr"^ + 21A^j^r + 2(86rj^A - 49A^ - 3hAyr)ry 
-AQvyyvAy + (85rj^A - 15 A, 
-((8A^ + 7 Ayr - 32ryA)r - 10r^A)r 
-((29A,. + 21Aj,r - 95ryA)r - 46r,A)r,.j,r - IQr^yyr^A 

+ 7^lil)yar - 6(^JaA + 7^yyyyr)<^yr^)<^l 
-(2(5((4A,, + 7Ayj,r' + 7A,yr - (23A, + 21Aj,r - 31r^A)7 

2 



— 17 ryyrA)r 
-10r,j,r2A) 



- ((9A, + 7Aj,r - 27rj,A)r - 6r,A)r, 
21(^j,V'j,yj,r5)(^J - 15((3((5A, + 7Aj,r 



- 3r^^r A 
llryA)r 



-4r^A) - 7ipyi)yyr^)ipy + 7 {if y^) yr " 6 A ) (/^yj^r ^ ) r ^ ) (/Jj^j^r 
-2((5((5A, + 7Aj,r - llryA)r - 4r,A) - U^yxljyyr'')^y 

+ b{7^yi)yr - 3QA)^yyr^)^yyy^yr^ + {7 ^ yi) yr " 3 A ) (/^yj^j^j, 



(39) 



{^lA)-' ((((2(r.., + 2r 



r")A 



+3(3A^ + 2Aj^r - 8ryA)ryyr^)r - ((lOr^. + llryr)A 
-(4A^ + Ayr)r)r^^ - ((13r^ + 20rj,r) A - (7A^ + 3Ayr)r)r, 
+ ((<^J/?V' + ^li'ya + ^j,yj,j;)r - ^laA)^yr^ + (9A^.^. + 4Aj^j,r2 
+7A,y 

~^ ^xxx 



2(13A, + 6Aj,r - l2ryA)ry)ryr^ - {{A,,y + Ayyyr^ 



+ A,j,yr2)r2)r - ((2((17A, + SA^r - 23ryA)ry + QryyrA) 



xxy 

yA)ry 



'■yyy' 
yy 



-(6A,, + Ayyr^ + 3A,yr)y - (5(3r, + 8r3,r)A 

-3(5A^ + Ayr)r)r^)r^)^l - ((2((5(r^^ + 3ryyr^ + 2r^yr)A 

+3Lpyiljyyyr^ + 5(5A^. + 3Aj^r - QryAYyr - 5(A^.^. + Ayyr"^ + A^yr)r)r 

-5((3r^ + 7rj,r)A - (3A^ + Ayr)r)r^)^l - 15((3(r^ + 2rj,r)A 

+^y'4)yyr^ - 3(A^ + Ayr)r)^y - {^y^yr - 7A)(fyyr^)(pyyr^)!fyy 

+ (2((5(r^ + 2rj,r)A + 2ipy^yyr^ - 5(A^ + Ayr)r)ipy 

-^{ipyipyr - QA)ipyyr'^)(pyyy ^ {^f y^^ yr " 5 A ) (/? y y f ^ ) ^9 r ^ ) r ^ . (40) 

Definition 2. We call fl2Tl) with arbitrary coefficients r = r{x,y), 
Fo = Fo{x,y), Fi = Fi{x,y), F2 = F2{x,y), Hq = Ho{x,y), Hi = Hi{x,y), 
H2 = H2{x,y), Jo = Mx,y), Ji = Ji{x,y), J2 = ^2(3;,?/), J3 = M^.y), J4 = Ji{x,y), 
and Ki = Ki{x,y), (i = 0, . . . , 7), the second candidate for linearization. 
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Thus, we showed that every hnearizable fourth-order equations belong either to 
the class of ([9]) or to the class of fl2T]) . In Sections 12.21 and 12.31 we formulate the 
main theorems containing necessary and sufficient conditions for linearization as well as 
the methods for constructing the linearizing point transformations for each candidate. 
Proofs of the main theorems and illustrative examples are provided in the subsequent 
sections. 

2.2. The linearization test for equation ^ 

Consider the first canditate for linearization, i.e. equation ([9]). In this case, the 
linearizing transformations ([2]) have the form 

t = V5(x), u = ip{x,y). (41) 

Theorem 1. Equation ([9]) 

^(4) ^A,y' + Ao)y"' + B^y"^ + {C^y'^ + C^y' + Co)y" 
+D,y'^ + Dsy'^ + D^y'^ + D,y' + = 0, 

is hnearizable if and only if its coefficients obey the following ten equations 

Aoy - Ai, = 0, (42) 

ABq - 3Ai = 0, (43) 

12^1^ + ^Al - 8C2 = 0, (44) 

12Ai^ + 3AoAi - ACi = 0, (45) 

32Coy + l2Ao,A, - I6C1, + 3AlA, - AAqCi = 0, (46) 

AC2y + - 24^4 = 0, (47) 

AC,y + A^C^ - 12D^ = 0, (48) 

16Ci^. - 12^0x^1 - 3AlAi + AAqCi + SAiCq - 32^2 = 0, (49) 

192D2^ + 36Ao^AoAi - ASAq^Ci - ASCq^Ai - 2SSDiy + ^AlAi - l2AlCi 

-36A0A1C0 + ASAqD2 + 32C0C1 = 0, (50) 



384Di,y 



3((3AoAi - AC^)Al + 16(2AiDi + C,C^) - IQ{A^C, - 1^2)^)^ 
-32(4(CiDi - 2^2/^0 + C0D2) + (3AiDo - Cl)A^) - 96Di,Ao 
+384Doy^i + 1536L>02,y - 16(3^0^1 - ACi)Cq^ 

+12((3AoAi - 4Ci)Ao - 4(AiCo - 4D2))Aox] = 0. (51) 
Provided that the conditions P2l) - fl5Tl) are satisfied, the linearizing transformation (j4Tl) 
is defined by a fourth-order ordinary differential equation for the function (p{x), namely 
by the Riccati equation 

40^ - 20x' = 8C0 - 2,Al - 12Ao,, (52) 
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for 

X = , (53) 

'^x 

and by the following integrable system of partial differential equations for ilj{x,y) 

ilPyy = iPyAl, (54) 

A^P^y = ^|JyiAo + 6x), (55) 

and 

1600^,,,, = 9600V'...X + 1607/-,, (-12^0. - 3^^ - 90x' + 8Co) 

+40iP,{12Ao,Ao + 72Ao.x - IGCox + SA^ + 18 Alx - UAqCq 
+120x^ - 48xCo + 24Di - SVl) + ^(144^^^ + 72AoxAl - 352Ao.Co 
-160^0^^ - SOCo^Ao - leOODoy + 640Di^. - 80fi^ + 9A^ - SSAICq 
+ mAoDi + 30 Aon - AOOAiDo + 300%^^ + 144Co2) + IGOOtpyDo, (56) 

where x is given by (!53|) and Q is the following expression 

n = Al- AAoCo + 8Di - 8Cox + GAo.Ao + AA^xx- (57) 

Finally, the coefficients a and (3 of the resulting linear equation ([3]) is given by 

and 

13 = {lQm^l)~\-lUAl^ - 12Ao^Al + 352Ao.Co + 160Co.. + SOCo.Aq 

+1600Dos/ - 640^1^ + - QAl + SSAlCo - imA^Di - ?,OAqVL 
+400AiL)o- 300x^^-144^2). (59) 

Remark 1. Since the system of equations fj42|) - fl5T]) provides the necessary and sufficient 
conditions for linearization, it is invariant with respect to the transformations ( HTl) . It 
means that the left-hand sides of (H2|) - (!5T|) are relative invariants (of the second-order) 
for the equivalence group (14T1) . 

2.3. The linearization test for equation (21\) 

The following theorem provides the test for linearization of the second candidate. The 
necessary and sufficient conditions comprise eighteen differential equations (!60 l) -( l77l) for 
twenty one coefficients of the ( 12T|) . The linearizing change of variables ([2]) is determined 
by (178|l - (l8T|) for the functions ip{x,y) and ifj^x^y). 

Theorem 2. Equation fl2Tl) 

y^'^ + FT7(-ioi/" + P^y'^ + Piy' + F^)y"' 



.- [15y"3 + {H^y''' + E^y' + Eo)y" 

H-hy'^ + -hy'^ + -hy'^ + Jiy' + Jo)y" 

+Ksy'' + K^y'^ + K,y' + K^] = 0, 
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is linearizable if and only if its coefficients obey the following equations 
IQryy = -{F^y + + Fa^^r + r^^Fa), 

lOr^; = IQryT - Fo + F^r - Fgr^, 

H2 — — 3F2, 

4/^1 = -3(5Fi -2F2r), 

4//o = -3(6Fo-Fir), 

lOFij,^ = -(Fij,F2 - 40F2,j, - 16F2,F2 + 2myyr + ^^yTy + UF2yF2r + 2QJ^ 

- 20 J^yV + UvyF^ - iOVyJ^), 

12^2^ = 12F2j,r - 3F1F2 + 6F|r + 4 J3 - IQJ^r, 

60Fi^ = GOFiyr - 36F0F2 - ISF^^ + 66FiF2r - 36F|r^ + 40 J2 - SOJ^r + SQJ^, 

60Fo^ = eOFo^^r - SIFqFi + 66FoF2r + 36FiV - 72FiF2r^ + 36F|r^ + 60 Ji 

- 80 J2r + 80 Jgr^ - 80 J^r^, 

20 Jo = 9F^ - ISFoFir + ISFoF^r^ + QF^r^ - 18FiF2r^ + 9F|r^ + 20Jir 

- 20J2r2 + 2QJ3r^ - 2Qjy, 

120 Jsyy = 216FiyF2y + ^AF^yF^ - 48Fi^ J4 + 360F2^j,r^ + 90F2yyFi - 180F2yyF2r 

- 432F|^r + 324F2j,r2;F2 + I89F2J/F1F2 - 486F2j,F|r - 192F2yJ3 
+ 86AF2yJ4r - 60J3yF2 + 120 J^y + 180J4^F2 - 240 J4j/yr 

- 1200 J42;r2; + 60 J4j;F2r + 720^60: - 120K^yr - 5040X7^r 

+ hOmK^yr"^ + 36rj,F2' - 432rj,F2 J4 - 2160ryi^6 + l^l20ryK-jr 
+ 504Foir7 + 36FiF| - IO2F1F2 J4 - hO\F^K^r - 12F^r 

- 48F2V3 + 396F2V4r + 50AF2K7r^ + 136 J3 J4 - 5447^, 

2A0J4xyy — ~{^^FiyF2yy + 162FiyF2yF2 — 72FiyJ4y + 36Fij^F| — I68F1J/F2 J4 

- 72FiyKe - 168FiyK7r - 72F2yyF2yr + UAF2yyryF2 

+ 54:F2yyFiF2 " 108F2j,j,F2V - 72F2yyJ3 + 288F2yyJ4r + A32F^yry 
+ IO8F4F1 - 540F4F2r - 144F2^ Ja^ + 528F2^ J4, + 192F2^ J4^r 
+ 324F2j,rj,F2' - 1008F22/r^ J4 + 162F2yFiF^ - I32F22/F1 J4 

- 396F2yF2V - 180F22,F2 Ja + 1320F2j,F2 J4r + 144F2^X6r 

- 336F2yir7r^ - 36Ja2;F| + 17673^^ J4 + 120 J4xj,F2 + 132J4^F| 

- 432J42; J4 — 240J4j^yj,r — 960J42^j^r2^ — \2OJiyyF2r — 768 J4yryF2 

- 138J4j;FiF2 + 288JiyF^r + 184J42;^3 - 1008 J4j;J4r + ^QOK^^y 
+ 240K6^F2 - 960Keyyr - SSAOKeyVy - 240KeyF2r - 1920K7^yr 
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+ 2AmQK^yryr + 432X72^^0 + l&^K^yFir + %l2KyyF2r'^ 

+ 20160rjir7 + 1728rj/FiX7 + SGr^/Fs^ - 264rj/F| J4 - l2A%ryF2KQ 

+ 5280rj;F2X7r + IGOr^, J| + 4O8F0F2X7 + 150^1^X7 + 21F^F^ 

- I2OF1F2V4 - I68F1F2X6 + 168FiF2X7r - 54F|r - 36F|J3 

+ 384F| J4r + 336^2^X6^ - 1344F2^X7r^ + I6OF2 J3 J4 - 64OF2 J4V 

- 400 J2X7 + 224 J3K6 - 368 J3ii:7r - 896 J4X6r + 3872 J4X7r2 
+ 672Fo,i^7), 

4J4^ = 4J42^r - Fi J4 + 2F2 J4r - 4X5 + 2AKQr - 84i^7r^ 

mF^yy = -{30FoyF2 + 36FiyFi - 36FiyF2r - 60F2yyr^ + 24F2j,Fo - 36F2yFir 

- bAF2yF2r^ - AQJ2y + 40 Jgj^r + 80 J4^r2 - 36ryF^F2 + 36ryF^r 
+ AOryJ^ - SOryJir + GFqF^ - 6F0 J4 + QF^Fs - 18FiF|r 

- 12Fi J3 + 24Fi J4r - 6F2V - IOF2 J2 + 22F2J3r + 26F2jy 

- 6QK4 + ISQK^r - ISQKer^ - mKrr^), 

20J2^ = 20J22,r + 20J3^r - 20J32,r2 - MFq J3 + 28Fo J4r - 5Fi J2 + IdF^J^r 

- 28Fi J4r^ + IOF2 Jar - 24F2 J3r2 + 28F2 J4r^ - I2OX3 + 360ir4r 

- 640X5r^ + 840^6^^ - 840 K^r^ 

60Ji^ = GOJiyV - iOJ^y + iOJsyr^ - 42Fo J2 + 42FoJ3r - 70FoJ4r2 - 15Fi Ji 
+ 42Fi J2r - 52FiJsr^ + lOFiJ^ + 3OF2 Jir - 42F2 Jar^ 
+ 62F2 J3r^ - 7OF2 J4r^ - 6OOX2 + 1080X3r - 1380X4r^ 
+ 1700ir5r^ - 2100^6/ + 2100i^7r^ 

8QK1 = 3F^Fi - 6F^F2r - 6FoFfr + 18FoFiF2r^ - 12FoF|r3 - SFq Ji 

+ IGFoJar - 24FoJ3r2 + 32FoJ4r3 + SF^^r^ - 12F^F2r^ + IhFiF^r^ 
+ 8F1 Jir - I6F1 J2r2 + 24Fi Jgr^ - 32Fi J4/ - 6F2V - 8F2 Jir^ 
+ I6F2 J2r^ - 24F2J3r^ + 32F2Jir^ + imK2r - 2Am^r'^ + 320i^4r' 

- 400is:5r^ + 480X6^-^ - ^mK^r^ , 

400X0 = -(6Fo^ - SSF^Fir + 48Fo^F2r^ + 48FoFiV^ - 126FoFiF2r=^ + 78FoF|r^ 
+ 40Fo Jir - 8OF0 J2r^ + 120Fo J3r2 - I6OF0 J4r^ - 21FiV=^ 
+ 78F^F2r^ - 93F1F2V - 40Fi Jir^ + 8OF1 J2r^ - 120Fi J3r^ 
+ I6OF1 J4r5 + 36F2V^ + 4OF2 Jir-'^ - 8OF2 J2r^ + I2OF2 J3r5 

- I6OF2 J4r'' - 4007^2^^ + 800K3r^ - 1200K4r^ + IGOOK^r^ 

- 2000 KeT^ + 2A00Ky). 
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Provided that the conditions (H7]) - fl77p are satisfied, the transformations ([2]) mapping 
equation fl2T]) to a hnear equation ([3]) is obtained by solving the following compatible 
system of equations for the functions ip{x,y) and ip{x,y) 

= r^py, (78) 
Vy-ipx = ripytpy - A, (79) 

lOAiPyy = ^y{AAy " ^2 A ) , (80) 

and 

bOOify^yyyyA^ = 300'^yyyiPyA'^{4Ay " FzA) + St/^j^j, A ( " 1 20^2 A^ - 144AJ 

+ 72AyF2A - 39F|A2 + 80J4A^) + ijyipy{-500iplaA^ 

- l50F2yyA^ + 360F2yAyA^ - 165F2yF2A^ + lOOJ^yA^ + 96Aj 

- I2AIF2A + 108Aj^F|A2 - 240Aj^ J4A2 - 24F|A3 + 6OF2 J4A3) 

- 500V'<^J/3A3 + 500i^7A^ (81) 

The coefficients a and (3 of the resulting linear equation ([3]) is given by 

« = A , (82) 



and 



(3 = (leOOAyp^)-^ [A(-144F4 - 72^22,^2^ + 352^2^74 + 160 J4yj; + 8074^,^2 

+640fs:6y - leOOfsTyx - 2SmKjyr + 80ej, - AASQryK^ - AmF^Kj 

-9F^ + 88F| J4 + l60F2Ke - 320^2/^7^ - lAAJ^) - 120Aj,e] , (83) 

where is the following expression 

= (F| - 4J4)F2 - 8(7^6 - 7Krr) - 8 J4, + 6^2,^2 + 4^2,,. (84) 

Remark 2. The equations fl60|) - fl77|) define eighteen relative invariants of the third- 
order for the general point transformation group ([2]). 



3. Proof of the linearization theorems 

The proof of the linearization theorems formulated above requires investigation 
of integrability conditions for the equations given in Section 12.11 We will 
consider the problem for the candidates ([9]) and ( 1211) separately. The problem is 
formulated as follows. Given the coefficients Ai{x,y), Bi{x,y),Ci{x,y), Di{x,y) and 
Fi{x,y), Hi{x,y), Ji{x,y), Ki{x,y) of the equations ([9]) and ( 12T1) . respectively, find the 
integrability conditions of the respective equations for the functions ip and ip. 
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3.1. Proof of Theorem 1 

Let us turn to the proof of Theorem 1 on hnearization of (j9]). Namely, given the 
coefficients Ai{x,y), Bi{x,y),Ci{x,y), Di{x,y) of ([9]), we have to find the necessary and 
sufficient conditions for integrabihty of the over-determined system (fT0l) - (!20l) for the 
unknown functions ip{x) and ilj{x,y). 

We first rewrite the expressions ffTOj) and ffTTl) for Ai and Aq in the following form 

^yy = (85) 

/ (Q^xx + ^xAq) 
^xy = (86) 

Comparing the mixed derivative {ipyy)x = {4'xy)y, one arrives at (j42l) 

Aoy = Alx- 

Then fll2l) . (IT^ and (I14p are written in the form 

3Ai - ABo = 0, 

3Al - 8C2 + UAiy = 0, 

and 

12^1^ + SAoAi - 4Ci = 0, 

respectively. So that one obtains fH51) . flH|) and fllSl) respectively. Furthermore, ffT51) for 

Co becomes 

{UAo^xVl - 60yL + ^vlAl - S^lCo) 
^^^^ = 40^; • 



Differentiation of flHTl) with respect to y yields 

12Ao^Ai + 32Coy - l6Ci^ + 3AlAi - AAqCi = 0. 

Thus one gets P6l) . Therefore f|T6l) . f|T7j) and f|T8|) can be written in the form of (j4 
l8|) and fH9|) . respectively. 

One can determine a from (HM . as the following 

AAoxx + QAo^Ao - + Al- AAoCo + 8D1 



Since = V9(x), we have = yields ([50 



- 



UAlCi - 36/loAiCo + A8A0D2 + 32^0^1" . 



From (!20!) one finds 



XX 



+ 36Ao,(/?^^,Ao - 48Co,</^^^, - 120(^L^, + 360¥pL</'x^xx 
- 240v9^^v5^V'xxx - 18v3xxV5^V'i^^o + 48v3xxV5xV'xCo + "iOiplPip 
+ Uifil^/j^^Al - ?>2^1^Ij^^Cq + Sv^^^/'^Aj^ - 2Q^IiP^AqCq 

+ m^lij^D^ - m^li,yD^. (89) 



Linearization of fourth- order ODEs 



14 



Forming the mixed derivative {ipxxxx)y = i'^xy)xxx one obtains 

- 1800Aox(fxxAo - UAo^Lf^Al + 32AoxLpxCo - ASOCo^^ipx 
+ 2mCox^xx + IGOODoy^x - SOO^x^Al + UOOif^^AoCo 

- 2mip,,D, - 39ip,A^o + 208^,AlCo - mip.AoD^ 

+ 400(^,^1^0 - U4VxC^] . (90) 
Since (p = ip{x), we have Py = yields 



n — 3 

J-^lxy — 



[(3AoAi - 4C,)Al + 16(2AiDi + CoCi) - 16(AiCo - /^2)Ao]^o 
-32[4(CiDi - 2C2D0 + CoZ^2) + (SAiDo - Ci)A,] 
-96DiyAo + 384:DoyAi + 1536Doyy - l6{3AoAi - 4:Ci)Cox 
+12[{3AoA^ - 4Ci)Ao - 4(AiCo - W2)]Aox • 

From (187|) one can rewrite the representation for Cq upon denoting x = ^ leads to 
and the representations for ipyy and ■i/'ajy in the equations (ISBj) and (15^ become (15^ 
and fl55|) . Rewriting the representation for a from fl88l) in the form 



where 



n = Aj* - 4AoCo + 8D1 - 8C0.. + GAoxAo + 4Ao 



and thus /3 of flQOl) becomes 

P = {1600ipt)-\-UAAl^ - 72Ao.Al + 352Ao.Co + WOCoxx + SOCo.Aq 
+ I6OOD02; - 640^1, + SOn^ - 9A^o + 88AlCo - WOAqDi - 30^0^ 
+ 4OOA1D0 - 300x^^ - 144C2). 

Finally, one obtains fl89|) in the form 

1600^,,,, = geOOV^xxxX + 160V^,,(-12Ao. - 3^2 - 90x' + 8C0) 
+ 40^,(12Ao.Ao + 72Ao,x - + SAl + ISAlx - l2AoCo 

+ 120x^ - 48xCo + 24Di - 8^]) + ^^(144^^^ + 72^0,^^ - 352/lo,Co 
- I6OC0,, - SOCq^Aq - leOODoy + 640Di, - SOn^ + 9A^ - 88AlCo 
+ leOAo/^i + SOAon - AOOAiDo + 300x^ + 144C2) + leOO^yDo- 
Hence we complete the proof of Theorem 1 . 

3.2. Proof of Theorem 2 

In the case of (pTj) . the problem is formulated as follows. Given the coefficients 
Fi{x,y), Hi{x,y), Ji{x,y), Ki{x,y) of (^I^, find the necessary and sufficient conditions 
for integrability of the over-determined system of equations ( l22|) - fH0|) for the unknown 
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functions ip{x,y) and ipi^x^y). Recall that, according to our notation, the following 
equations hold 

^x = r^Py, ipx = , (91) 

and 

_fx R -!£^R 

(^x ^yy Px Py 

fy 'Py 
Let us simplify the expression fl2^ as follow 



'Pyy 



{AAy-F2A)^y /(lOA). (92) 



Comparing the mixed derivative {fx)yy = {fyy)x one obtains 



^xy 



F2,A^ - F^yvA^ + lOryyA^ + 4ryAyA - r^^FsA^ 

.2 ' 



+ 4A,Aj, + 4A^j,rA - 4A^rJ /(4A). (93) 
Rewriting (123|) in the form 

A^ = (20rj^A + 4Aj^r + FiA - 2F2rA)/4. 
Forming the mixed derivative A^y = (A^)j^ one arrives at (1601) 

'^yy = ~{Fiy - F2x - F^yT - rj/F2)/10. 

Then f l2^ - fl27|) are written in the form of fl6Tl) - flMl) . respectively. Furthermore, fl28l) 
becomes 

Ayy = -{20F2yA^ - A8Al + AAyF^A + 7F^A^ - 20JiA^)/{A0A). 

Now, consider the equation {Ayy)^ = {Ax)yy , one gets (l65l) 

F^yy = -{F,yF2 " 40^2,^ " 16^2,^2 + 20F2^j,r + 40F2yry 

+UF2yF2r + 20 - 20 J^yr + UryF^ - 40rj^ J4)/10. 

Thus equations ([29])- ([32]) yield ([66]) -([69]), and from ([33]) one finds 

i^yyyy = [SOOij yyy^ y A^iAy " ^2 A ) + 5^^j, A ( " 1 2 0^2^ A ^ - 144A2 

+ 12AyF2A - 39F|A2 + SOJ^A^) + i)y<^y{-hm<^laA^ 

- l^0F2yyA'^ + 360F2yA^A2 - 165F2j,F2A3 + 100 J4yA3 

+ 96Aj - 72AJF2A + 108Aj,F|A2 - 240Aj, J4A2 - 24F|A=^ 

+ 6OF2 J4A^) - 500^(^J/3A^ + SOOiTrA^] /(SOO^^^A^). (94) 

One can determine a from fl34]) . as the following 

a = {4F2yy + 6F2yF2 -8J^y + F| - 4F2J4 - 87^6 + 56Krr)/8y^l. (95) 

Now the equation — ray = leads to fITDl) . Furthermore, one considers {ipx)yyyy = 

i'^yyyy)x , yields 

f3 = 120Ay{-AF2yy - 6F2yF2 + 8J4y - F| + 4F2 J4 + SKq - b^^r) 

+A{?>20F2yyy + 480F2yj,F2 + 336F|^ + imF2yFl + 32F2y J4 .gg^ 
-480 J4TO - 240J4yF2 - 16007^7^ + 1600ir7j,r - AmF^Kj ^ ' 

-9Fi + 88FIJ4 + 160^27^6 - 320F2K7r - 144J|)/1600A^^. 
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The equation - r/5j, = leads to Therefore, dSSD-dinD become ([72]) -([77]), 

respectively. 

Let us turn now to the integrability problem. One can find all fourth-order 
derivatives of the functions and ip by using (19T]) . ([92]) and ([Mj) . So that one obtains 
at ( |78l) - (!8Tl) . Finally, the coefficients a and (3 of the resulting linear equations ([95]) and 
f [96]) are given by 

e 



/? = (-144F4A - 72F2yF|A + 352^2^, J4A + 160 J4toA + 80J4j,F2A 
+640fs:6|/A - 1600is:7^.A - 2SmKjyr^ - USOryKjA + SOOyA 
-120Aye - 4OOF1K7A - 9F2^A + 88F| J4A + I6OF2KQA 
-320F2K7rA - 144J|A)/(1600(/?^A), 



where 



e = (F| - 4J4)F2 - 8(7^6 - 7Kjr) - 8J4, + 6^2,^2 + 4^2,,. 
Hence we complete the proof of Theorem 2. 

4. Illustration of the linearization theorems 
4.I. An example on Theorem 1 

Example 1. Consider the nonlinear ordinary differential equation 

x^y{2y^^'^+y)+Sx%'y"'+lQxyy"'+Qx^y"'^+ASxy'y"+2Ayy"+2Ay'^ = 0.(97) 
It is an equation of the form with the coefficients 

Ai = — , Aq = — , Bq = — , C2 = , Ci = — , Cq = — , 
y X y xy x^ 



D, = 0, Do = l 
x'^y 2 

One can check that the coefficients ([98]) obey the conditions ([^2]) - ([5T]) . Thus, the 
equation ([97]) is linearizable. We have 

8C0 - 3Al - 12Aq^ = (99) 

and the equation ([52]) is written as 

ax 

Let us take its simplest solution x = . Then invoking ([53]) . we let 

ip = X. 

Now the equations ([54]) - ([55]) are written 

Ipyy 1 Ipxy 2 

2/' X 
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and yield 

Tpy = Kx'^y, K = const. 

Hence 

i^ = K^ + fix). 
Since one can use any particular solution, we set K = 2, / (x) = and take 

; 2 2 

ip = X y . 

Invoking fl99l) and noting that fl571) yields = 0, one can readily verify that the function 
ip = x^y^ solves equation ( !56l) as well. Hence, one obtains the following transformations 

t = x, u = x^y"^. (100) 

Since i7 = 0, equations (158|) and (!59|) give 

a = 0, /5 = ^ = 1 

Hence, the equation (!97|) is mapped by the transformations fllOOp to the linear equation 

U^^) + M = 0. 

Example 2. The third-order member of the Riccati Hierarchy is given by Euler et al. 
[5] as 

y'" + Ayy" + ?,y'^ + Gy^^/' + Ay^ = 0. (101) 

Applying [6], and [7] one checks that equation cannot be linearized by a point 
transformation or contact transformation or generalized Sundman transformation. 
Under the Riccati transformation y = — the equation f llOip becomes [8J 

cu=^cu(^) + 4 (a - 1) cu^cu'cu'" + 3 (a - 1) uj^uj"^ 

+ 6 (a - 1) (a - 2) cjcj'V + (a - 1) (a - 2) (a - 3) cj'^ = 0. (102) 
It is an equation of the form ([9]) with the coefficients 

^ = ^, ^0 = 0, So = ^, 

^ ai^6a!+lla^ ^ 1^3 = , ^2 = , Z^l = , = 0. 

One can verify that the coefficients (11031) obey the linearization conditions (H2]) - fl5ip . 
Furthermore, 

8Co - 3Al - 12Ao^ = (104) 
and the equation ( l52l) is written as 

ax 

We take its simplest solution x = and obtain from (153|) the equation ip" = 0, whence 

If = X. 
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Equations and have the form 



and yield 



Hence 



^/^^ = fsTcj^'^-^), = const. 



iP = K— + f{x). 
a 

Since one can use any particular solution, we set K = a, f (x) = and take 
^ = uo". 

Invoking (11041) and noting that fl57|) yields = , one can readily verify that the function 
ip = uj"" solves equation flS^ as well. So that one obtains the following transformations 

t = x, u = uj''. (105) 

Since f2 = 0, equations (158|) and (!59|) gives 

a = 0, /3 = 0. 

Hence, the equation (I102p is mapped by the transformations (llOSp to the linear equation 
w(^) = 0. 

Example 3. Let us consider the Boussinesq equation 

utt + uu^^ + ul + u^^^^ = 0. (106) 

Of particular interest among the solutions of the Boussinesq equation are travelling wave 
solutions: 

u{x,t) = H{x - Dt). 
Substituting the representation of a solution into (]106p . one finds 

H^*^ + {H + D^)H" + H'^ = 0. (107) 

It is an equation of the form (Q with the coefficients 

A, = 0, Ao = 0, Bo = 0, C2 = 0, C, = 0, Co = D' + H, 

^4 = , 1^3 = , 1^2 = 1 , ^1 = , Do = 0. ^ ' 

Since the coefficients fllOSp do not satisfy the linearization conditions ( H6ll . (H9l) and 
flSTj) . hence, the equation (11070 is not linearizable. 

Example 4. Consider the non-linear equation 



^y"y"' + \ (15/^ - xy'' - y") = 0. (109) 



y u 
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It has the form (PT]) with the following coefficients: 



r = , F2 = , 


i^i = 0, 


Fo 


= 0, 


H2 


= 0, H, 


= 0, Ho 


= 0, 


J4 = , J3 = ^ 


, ^2 = 


, ^1 


= 0, 


Jo 


= 0, Kr 


= -X, 




i^6 = -l, i^5 = 


= 0, 


= 0, 


. K3 - 


= 


, K2 = 




, Ko = 



Let us test the equation fllOQp for linearization by using Theorem 2. It is manifest that 
the equations ( l60l) - (177l) are satisfied by the coefficients f lllOp . Thus, the equation (11091) 
is linearizable, and we can proceed further. 

Let us take its simplest solution ip = y and ip = x which satisfy the compatible 
system of equations fl78|) - fl8T]) . So that one obtains the following transformations 

t = y, u = X. (Ill) 
Since 6 = 8, equations ( 182|) and ( l83l) give 
a = l, p = l. 

Hence, the equation fll09p is mapped by the transformations fillip to the linear equation 

M^^) + U + U = 0. 
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